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Contest 4 


In the diagram on the cover, the consecutive 
integers are entered into a spiral of squares, 
starting with one at the center, and proceeding 
clockwise. 

For every prime number (circled) entered, a 
square is skipped. If a prime falls next to a 
previous prime (as happens four times up to 100), 
two squares are skipped. (The "next to" is taken 
only horizontally or vertically, not diagonally.) 

The sequence 1, 5, 14, 30, 54, 84,... (that is, 


the numbers in the squares that proceed northeast 


$25, or a two-year subscription. 


from the center square) is to be extended, The 
first (and only) prize in this, our 4th contest, 


goes to the one who produces the longest correct 


The prize: 


list of the contents of those squares, done by 


computer. All entries must be received by POPULAR 


COMPUTING by April 30, 1976. (i 


POPULAR COMPUTING is published monthly at Box 272, Calabasas, California 91302. Subscription 
rate in the United States is $18 per year, or $15 if remittance accompanies the order. For Canada and 
Mexico, add 54 per year to the above rates. For all other countries, add $6 per year to the above rates. 
Back issues $2 each. Copyright 1976 by POPULAR COMPUTING. 


Publisher: Fred Gruenberger Contributing editors: Richard Andree Advertising Manager: Ken W. Sims 
Editor: Audrey Gruenberger William C. McGee Art Director: John G. Scott 
Associate Editors: David Babcock Thomas R. Parkin Business Manager: Ben Moore = 


Irwin Greenwald 
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The function defined as follows: 


PC35-3 


The CSR Function 


(Continued Square Root) 


A= 


(see PC34-10) provides some interesting computations, 


Q, + GQ, + GQ, + O, AE Je a | 


using nothing more advanced than square root and addition. 
The accompanying table shows the value for A for various 
sequences of a's, calculated by Herman P. Robinson on 

his Wang 720C. 


Consecutive integers ls 


Den integers 


Odd integers 
Fibonacci 
Powers of 2 
Squares 
Cubes 

4th powers 
5th powers 
10th powers 


all ones 


1,4, 2, 8, 5,7 repeating 


Primes 


Aaron 


2 
at 
aL 
al 
1 
2 
2c 
a 
6 
1 
1 
2 


75793 27566 


.15847 68723 
.85025 31288 
66198 24623 
-78316 58092 
94265 54227 
-17678 85971 


46740 45317 


.82348 15128 
.05181 33295 
-61803 39887 
-87349 51093 
- 10359 74963 


Herschfeld, 


Monthly, Vol. 42, 1935, 
sequence will converge if 


18004 
11039 
25914 
27811 
64098 
63987 
33441 
Lass 
34203 
21526 
49894 
71315 
39897 


53270 
76565 
28891 
55796 
88271 
32822 
98583 
10674 
37757 
60210 
84820 
48791 
26261 


88196 
58534 
21451 
76060 
03049 
14132 
98570 
79906 
13444 
47308 
45868 
93475 
99396 


38218 
79807 
aio 
81815 
92255 
91412 
61721 
49563 
04555 
44812 
34365 
30993 
49685 


13852 
02524 
Oa 
13129 
00328 
66723 
83836 
77992 
84089 
93383 
63811 
64755 
32544 


in the American Mathematical 
pp. 419-429, showed that a 


lim (1n In PEON Aeabel )ikge ees) 


n-*o 


a,+ ¥a4-- 


76531 
16696 
79416 
50561 
85836 
76880 
53483 
98547 
68438 
THING 
77203 
34321 
4o4e21 
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Zhe 
Thus, a sequence with terms increasing no faster than xe 
will converge. al 

If the terms are exactly x“ , the CSR is x7, where 7 is 
the golden mean. The sequence is started with n=l. 


Mr. Robinson points out the paradoxical situation 
that as the terms grow faster in size, fewer terms are 
needed for convergence, in general, An example illustrates 
this: 


Let a given term be Tose. The following term will 


be Patel”. but the square root will be Vk"10°, Dropping 


this term introduces an error of the order of Vk parts in 
io?” if k is not too large. Successive square roots 
further reduce the error drastically. In the case of the 


factorials, starting with n = 17 gives an answer with a 
fractional error of the order of BmGr ad Robinson 
estimates that starting the factorial calculation with 35! 
will give a result good to more than 290 decimal places. C] 


Log 35 1.544068044 3502756 35498477 36386814 3166715382514861857 


nn 
© 
U) 
LU V35 
ag 
WW 
0 
Z 


. 55534806 14894 1 3679706112076669 367 36916268608 3850379 
-91607978 3099616042567 328291561617 048415501230794 340 


3 

Wath AS 

Y35 3.271066 3101885897 282248069023925 31344098903147778906 

¥/35 2.036168004640398017 360874 164145317694261816167578535 
1.42694 358845765098 359004986 1650304 28887 1125929977205 
1.03619 306288839615 3570156125080192301740086857 188821 


e2? 1586013452313430 . 7281296446257 74660125176203950134526 
154266697022452801204626923251641062 


T 251330702007 364298 .616088947097 500656998 3291245887646 
0025 367 2035657 158609090814 25139385 


tan ~ 35 1.542232668956136624759150722706513354 393176599772952 


Many computing problems can be solved by means of 
e a method called bracketing. The method (like bubble 
sorting) is inelegant, crude, unsophisticated, and 
inefficient. On the other hand, it is delightfully 
simple, easy to code in any language, and applies to a 
broad range of problem situations. Moreover, it is 
widely used, and knowledge of it should be made available. 


It is best explained with an example. Suppose we 
are using a computer with a 32-bit word size. The high 
order bit of a data word is used to indicate the sign, 
leaving 31 bits for the size of the number. Thus, on 
such a machine, the largest positive integer that can be 
contained in one word is 2,147,483,647--one less than the 
31st power of 2. 


We wish to construct a subroutine that will accept 
a number, N, and produce in a word addressed at SR the 
square root of N. The range on SR is zero to 46340. 
A perfectly logical algorithm, then, would be to try 
every value between O and 46340 to find which one, when 
squared, came closest to N. If N is around two billion, 
the algorithm will require 46000 traverses of its loop in 
the worst case, which is hopelessly inefficient. 


We should take bigger steps to start with, and then 


cut down the step size as we approach the root. The 
scheme shown in flowchart A will do nicely, keeping in 
mind that we are dealing entirely in integers. The 


value picked for the initial step size, 32768, is the 
largest possible power of 2 within the permissable range 

of SR. Since we choose to cut the step size in half 

each time we pass Reference eis we will pass through PE 
Reference 3 not more than 15 times to locate any root 

to the nearest integer, 


This is the simplest application of the bracketing 
process (and there are, of course, much better algorithms 
for square rooting). We can see the essential elements 
of the process, however, from the simple example: 


, e 
eo ART «- COMPUTING 10 —=s 


*The method resembles interval-halving here, but the two 
should not be confused. -In general, intervali-halving 
procedures are more difficult to program. A discussion 
of interval-halving can be found in Introduction to 


@ Computer Science, by Harry Katzan, Jr., Petrocelli/Charter 
books, 1975. 
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1. The solution space should be known; that is, the 
total range within which the result must lie. 


2. We arrange to explore that space, starting with 
very large steps. 


3. The key element is the decision (Reference 4 in 
flowchart A). We must have a way to determine when we 
have gone passed the desired result. Notice that this 
matter is not self-correcting. If the result is passed 
and that fact is not correctly noted, the process will 
hang up and/or yield erroneous results. When the decision 
says that we have gone too far, we back off one full step. 


4, We arrange to cut down the step size by any 
convenient amount. (Most frequently the step size is 
reduced by a factor of 10.) 


5. We arrange a graceful way to terminate the 
process when we have reached the desired level of precision. 


Like all numeric processes, bracketing can not be 
applied blindly. The details will differ from problem 
to problem, but the overall plan is as stated. 


We now have a new tool, and we have tried it out on 
a known case. 


Suppose now that we are working in scientific 99 
notation, so that the range on N is from zero to 9:10°~. 
We could apply the same procedure, but our initial step 


size would have to be something like 1099, and in the 
worst case we would go around the loop 332 times. The 


bracketing process does not apply in this case. 


We can, however, apply it to other problems quite 
profitably. To find the real root of Wallis' equation 


{x2 - 2X - 5 = 0) for example, we have 
Ne =X) ON 5 


and the flowchart logic B applies. Preliminary analysis 
tells us that the real root lies between 2 and 3, which is 
our solution space. A counter, K, will terminate the 
process after 8 traverses of Reference 3, for a solution 
correct to 8 significant digits. The termination 
procedure could be any test for convergence, but it is 
simpler to count through the loop. 
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The process lends itself to problems in maxima 
and minima. For example, consider the old problem of 
constructing, from a sheet of tin of length L and width 
W, an open topped box of maximum volume. The logic of 
flowchart C shows the bracketing process applied to this 
familiar calculus problem. The volume calculated for 
each trial value of X is compared to the previous volume, 
gv. As long as the volume keeps increasing, we are 
heading in the right direction. As before, when the 
volume drops ("too far''), we back off one full step. 
But now there is a new problem. The various cases that 
exist in seeking a maximum or minimum can cause trouble. 


For safety, the box at Reference 3 should be made 
X - 2D—~- xX, to back off two full steps. 


Figures P, Q, and R show what can happen in seeking a 
maximum (indicated by the arrows). P is the normal 
case, in which a drop in the calculated value indicates 
correctly that the maximum has been passed; in this 
case, a retreat of one step would suffice. But in the 
situations shown at Q and R, although the drop does 
indicate that the peak has been passed, a retreat of 
one step would not bring us past the peak, and the 
process would be out of control. 


Again, we have applied the new tool to a problem 
that is better handled with calculus. But the new tool 
now allows us to operate easily in cases where a solution 
by calculus might be awkward. Consider, for example, 
the problem given in issue No. 16, page 15: To divide 
the number 10 into two parts so that when the first part 
plus its square root is multiplied by the second part 
pius its square root, the product 1s 23; that is, to solve: 


(xt JR (Go. x Vior= x =e, 


There are efficient schemes for solving such equations, 
but bracketing will also work, and it provides a quick: 
and dirty solution. 


A similar situation is found in the Lake/Fence 
problem (issue No. 15, page 10), where it is required to 
find the maximum of 


A = x(500 - 2X) - 5280°T + H(250 ~ x) 


where T = tan“+}(250 - x) /i » 


and H°= 27815900 + 500X - x* 


PC35-9 


*CuUTXeu UF weTqorid atdwts e *‘xoq uT3 peddoq-usdo 
oud JO wetTqoad 9y9 09 paettTdde ssssoad Bupqyayoeaq sau, 


Xia 


° 


(cae M) (xe Sore 2 AN 


” 
“ 


- 


rad 


A founToa 
ageTnopteg 


(ae AG <—A 


OT-SE0d 


The functional values 
being calculated and 
compared. 
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The bracketing process can even be applied (although 
it is tricky) to more complicated situations involving 
several variables. The problem of Snoopy's water dish, 
for example (see Figure S) calls for the dimensions of 
a dish in which the length of the cross section (E-F-G-H 
in the Figure) to be 24 inches, and the length X and the 
angle A are to be determined so that the dish will have 
maximum volume. 


P Blips 2X fe} 


The bracketing process is most useful for inverting 
a function. To take a simple case, suppose we need the 
arccosine function and we are working in Fortran, which 
gives us the direct function, cosine, but not the inverse. 
Flowchart D shows the logic of a subroutine for which the 
input is a number A between + 1, and the output, xX, is 
the angle whose cosine is A. (Flowchart D is modified 
from one made by Robert White). 


The subroutine logic begins by editing the input 
data for validity and then immediately takes care of the 
obvious cases. The bracketing process begins at Reference 
2, and the critical element is the comparison of cosine X 
with A at Reference 5. The test at Reference 4 is 
inserted as a safety factor, to insure that the process 
does not exceed the range of the principal value for 
arccosine. The scheme shown could be further improved 
by capitalizing on the symmetry of the cosine function 
in the first two quadrants. It might be further improved 
by testing for other common values, such as 0.7071067 
or 0.5 for A, giving "exact" values of 7/4 and 17/3 
respectively. 


The function arccosine can be calculated more 
directly by formula, of course, but the algorithm applies 
to any function. For a function that is analytic, there 
are usually analytic methods available to invert it, but 
the bracketing process can invert a function that is 
discontinuous. One can invert a life insurance premium 
table, or a postal rate scale, or almost any function 
that is defined in one direction. faa] 
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The game of Master Mind has been featured in 
Games & Puzzles since 1973, and was written up in Time 
in the December 1 issue. 


Master Mind is related to the game of Guessword 
Puzzles, which first appeared in Computing News in 
October, 1954. Guessword Puzzles is a game for two 
people, played with pencil and paper (a convenient grid 
for the game is given on the facing page). 


One player selects a (hidden) 7-letter word. The 
other player then guesses any number of 7-letter words, 


one at a time. For each word he guesses, the first 
player tells him how many letters agree with those of the 
hidden word in the same position. From the numbers that 


the second player receives from the first, he must 
deduce the hidden word. 


For example, suppose the hidden word is SUCCUMB. 
If the second player guesses successively GENERAL, 
ELEMENT, ADAMANT, and BANANAS, he will be given counts 
of zero in each case. The situation is now: 


aan 
LAL DIA [m4 


Guessword Puzzles 


The spaces at the top of the grid provide room to record, 

for the second player, the letters he has so far eliminated. 

If he now guesses UNUSUAL, he gets a count of one. Since 

the AL in UNUSUAL has been eliminated by the AL of GENERAL, 

he Knows that he has located a correct letter within the 

first five letters of the hidden word. A guess of UNTRIED 

(with a count of zero) further narrows the choice, and a 

guess of TRUSTED then isolates the second U of the hidden > 
word. The situation is now: 
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At this point, most of the available 
information has been used up, and 
the game continues. If the next 
guess was ICICLES (for a count of 
one) a good player would probably 
deduce SUCCUMB in four more guesses, 
for a score of l2 for the game. 


Since the game is active for one player and passive @ 
for the other, it becomes expedient to play two games 
simultaneously, and the playing grid is constructed for 
that purpose. Par for the game is about 20 guesses, 
and a game usually takes about 45 minutes. 


The two lists below are guesses for two different 
games, which refer to two different hidden words. The 
count for every word in the list on the left is one; the 
count for every word in the list on the right is zero. 
What are the two hidden words? 


ICICLES MINIMUM 
MINIMUM ICICLES 
SUPPORT GENERAL 
GENERAL ELEMENT 
UNUSUAL UNUSUAL 
ADAMANT FUTURES 
UNIFORM ODOROUS 
DIAPERS ADAMANT 
CRANIUM BOLOGNA 
PRESORT BANANAS 
BANANAS QUICKLY 
AVERAGE _ OBLOQUY 
FUTURES DIAPERS 


ABS OLVE ROADMAP CL] 
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= Problem: A subroutine in a program produces as its 
Output 3-digit numbers (N) in no special order and 
with possible strings of duplicates. A second 


subroutine is needed which will examine these 
numbers as they are produced and report each new 
longer string of duplicates. 


The flowchart shown below gives the required logic. 

A counter, C, keeps track of the string lengths. The 
counter is initialized to zero in the housekeeping 
phase of the main routine, and the values of old N (@N) 
and old C (ZC) are also initialized to zero. The 
second subroutine is called each time a new value of N 
is generated by the first subroutine. 


To test this logic, generate the sequence 


ipa a Sir 36 3,4,4,4,4,5,5,5,5,5,4,4,4,6,6,6,6,6,6, oO 
and verify the following expected output: 


ON Fw 
AN FWh eH 


and so on L] 


Enter : 1 Return 


Counting String Lengths 
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